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I. INTRODUCTION 



In recent years, the covariant approach to Lorentzian loop quantum gravity has converged 



on a single model - the "new" spinfoams [l, 2|. The model is crucially based on the embedding 
of unitary SU{2) representations into unitary SL(2, C) representations. In retrospect, the 
use of unitary SL(2, C) reps seems unavoidable: just as the quantum mechanics of particles 
naturally leads to unitary representations of Poincare, so should the quantum mechanics of 
spacetime elements naturally lead to unitary representations of Lorentz. The unitary irreps 
of SX(2,C) in the principal series are labeled by {p,j), where p is a real number, and j 
is a spin (a non- negative integer or half- integer) . The evidence strongly suggests that the 
representations which should be used are of the form j), where the Immirzi parameter 7 
is an arbitrary real constant. Heuristically, 7 can be said to arise from an axial term which 
can be added to the Einstein-Hilbert action without changing the classical field equations: 

S Ho ist[e, 00} = J ((e A e)* + 7~ 1 (e A e)) A F[u\ . (1) 

7 is conventionally taken to be positive. Nevertheless, negative values are just as legitimate, 
and we will consider them below. It appears that the only special values of 7 are and ±00, 



the latter corresponding to the Barrett-Crane model J3]. Furthermore, it appears that a 
fully satisfactory construction is only possible for finite 7. In particular, the Barrett-Crane 
model does not capture the angular degrees of freedom in the graviton propagator j4|. 

We wish to study the parity transformation properties of the Lorentzian model with 
finite 7. Since there are no signs of trouble with (C)PT invariance, in what follows we will 
consider P- invariance and T- invariance as equivalent. At least naively, the quantum theory 
with a fixed finite 7 is not P-invariant. This can be seen from ([T]), or from the fact that the 
(p, j) representation with p 7^ is chiral: it relates e.g. tx boosts with yz rotations in a way 
which requires the right-hand rule. From the point of view of the action ([1]), the chirality 
of the model is a quantum effect. One can therefore hope that in the large-distance limit 
corresponding to semiclassical GR, parity will be restored. This argument, however, is far 
from robust: 7 enters crucially into the construction of spacetime itself, rather than just 
providing "quantum corrections" to its dynamics. 

Another naive expectation from the action (CQ) is that a parity transformation should 
be equivalent to the replacement 7 — > —7. In the quantum theory, this statement must be 
handled with care. 7, and with it potential parity violation, enters the theory in two separate 



places. First, it enters in the vertex amplitude A v (hi) associated with the SU(2) elements 
hi on the links of a boundary graph. Second, it enters through the canonical commutation 
relations into the interpretation of boundary states ij){hi) as 3d geometries. As we will see, 
it appears that the second role is the one leading to problems with parity. 

A good acid test for the P-invariance of the theory at large distances is the calculation 
of the graviton propagator js-8]. The current spinfoam model was motivated in part by 
the desire to get a spin-2 graviton [9j. This seems to have been achieved , llOfl . and 
it is now time to worry about the graviton's chirality. A difference between left-handed 
and right-handed gravitons, like any deviation from GR at large distances, would be highly 
problematic for two reasons. First, we of course wish to reproduce classical GR with its 
empirical success. Second, on the interacting level, GR is the only unitary and Poincare- 
invariant theory for low-energy elementary spin-2 particles. Therefore, once gravitons are 
present, the agreement of their behavior with GR expectations is a necessary condition 
for the very existence of a consistent flat-space limit. Note that neither of these issues 
can be settled at the propagator level: the propagator can always be made non-chiral by 
rescaling the graviton field components. However, this would constitute a modification to 
the geometric interpretation of boundary states. 

Some information on parity issues is already available for the Euclidean version of the 
theory. First, the Euclidean EPRL vertex is manifestly invariant under the replacement 
7 — > —7. To see this, consider the vertex amplitude for a semicoherent state: 

J 1 i=± 

where j t = (|1 ± r y\/2)ji. Setting 7 — > —7 interchanges the and j\ factors, without 
changing the result. Also, the large-spin 4-simplex transition amplitude for a semicoherent 
boundary state is composed of two parity-related critical points, with equal weight to 
each (plus two unwanted terms which are peculiar to the Euclidean). 



On the other 



4-simplex limit [7], and found to violate parity: it is not invariant under odd permutations of 
the 5 nodes, which correspond to parity-odd 4d isometries. As could be expected, invariance 
is restored in the 7 = limit. The invariance of the vertex under 7 — > —7 implies that the 
problem arises from the geometric interpretation of the boundary state. It is conceivable 
that this chirality is an artifact of the 4-simplex approximation, and will become negligible 



rand, the Euclidean graviton propagator was calculated in the large-spin 



when larger graphs and spinfoams are taken into account. Still, the evidence is troubling, 
and it would be desirable to resolve the problem already at the 4-simplex level. 

In parity-related questions, extrapolation from the Euclidean theory to the Lorentzian is 
clearly problematic. We must address the issues in the Lorentzian case directly. Important 
analytical progress has already been made in this direction. The semicoherent transition 
amplitudes in the large-spin 4-simplex limit have been studied in [12I . Using the framework 
developed there, the graviton propagator in the same limit was recently analyzed 8j. Un- 
fortunately, with current analytical methods the Lorentzian theory is less transparent than 
the Euclidean. As a result, parity invariance has not yet been successfully addressed. This 
is the task of the present paper. 

In a nutshell, the parity features we find for the Lorentzian theory are analogous to those 
found in the Euclidean. In particular, we calculate the "semicoherent" piece of the graviton 
propagator, i.e. the piece that arises from variations of the rotation and spinor variables at 
fixed spins. This is the part of the propagator which is potentially problematic under parity. 
We calculate this piece for two parity-related components of the propagator, and find that 
the values differ by the same ratio e 2m ^ 3 as in the Euclidean 3]. Thus, the chirality problem 
with the graviton propagator persists in the Lorentzian model. 

Our study of the Lorentzian 4-simplex was done numerically, using a Python script. The 
script is based on the definitions in [12[, together with the formulas for the metric insertions 



from 



8] (there is a slight 



discrepancy in conventions between the two papers, and we stuck 



to the conventions of 12]). Some technical comments regarding the script are given in the 
Appendix. The script files themselves are included in the arXiv submission. 

In section [Til we demonstrate analytically that the Lorentzian vertex, like the Euclidean, 
is invariant under 7 — > —7. In the process, we find a phase discrepancy between the vertex 
as defined in terms of group characters (e.g. (l3|) and the one used in [12j. We argue that 



the version in 



13] is the more appropriate one. In section IIII1 we present numerical results 



for the 4-simplex. These include the chirality of the graviton propagator, as well as some 
symmetry relations between the parity-related critical points of the transition amplitude. 
These symmetry relations have a very simple form, which can certainly be derived ana- 
lytically. In section IHI C\ we present a minor correction to the geometric construction of 



semiclassical boundary states in 



12| . In section HVt we discuss future prospects. 



4 



II. THE 7^-7 SYMMETRY OF THE VERTEX 



A. Derivation from the group-character definition of the vertex 

The vertex amplitude A v as a function of SU (2) elements h\ on its surrounding graph is 
given by jl3|: 

AvM = [ TT d 9n IT E( 2 -?' + x ) 2 / dh x^ihgsmgm) , (3) 

where we made the dependence on 7 explicit in the notation, n labels the graph's nodes, 
I labels the links, and s(l),t(l) are the source and target nodes of the link I. As usual, for 
a graph with N nodes, we integrate over just N — 1 Lorentz group elements g n , with the 
remaining one fixed to the identity in order to make the integral finite. In this subsection, 
we will show that A vn = A„ _ 7 . For this purpose, it's convenient to think of the integration 
variables in fl3]) as actual SL(2, C) (or SU(2)) matrices, rather than abstract group elements. 
To streamline notations, we introduce a symbol for the J-conjugate of an SL(2, C) matrix: 

g J = JgJ- 1 = (g-y , (4) 

where the J in the second expression is the standard parity operation on 2-spinors. The 
operation (TJJ preserves the multiplication order and the Haar measure, and leaves SU(2) 
matrices unchanged. It can be interpreted as either a P or a T reflection of the Lorentz 
rotations. 

The inequivalent SL(2, C) matrices are fully characterized by their complex trace. In 
its eigenframe, such a matrix is seen to consist of a spatial rotation and a boost in the 
perpendicular plane. The trace then encodes the rotation angle and the boost parameter. 
The complex conjugate of a given trace encodes (for instance) the same rotation with the 
opposite boost. Intuitively, this corresponds to a symmetry of the unitary Lorentz reps: 

ti B = \FA x P ' J i B ) = X~ p ' j i A ) ■ (5) 

We will use below a special case of this relation: 

x p 'V) = x~ p ' 3 (g) ■ (6) 

Let us present a short proof of ([6]). The equivalence of matrices with equal traces will then 
imply the general relation ([5]). The (p,j) representation of SL(2, C) is given by functions of 
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a 2-spinor z with the homogeneity property: 

f(\z) = \-l+3+iP\-l-j+iPf( z ) } (7) 

with the group action and the Hilbert product defined by: 

U 9 [f](z) = f(g T z) (8) 
(A,/ 2 >= / SlUz)h{z) . (9) 

J CP 1 

The integral is over non-collinear spinors z, and the two-form measure Q is defined as: 

i 

Q = -(zodzi — zidzo) A (zodzx — z~\dz§) , (10) 

where Zq and z\ are the two components of z. Note that Q is real. Now, consider the 
representation obtained by acting with the matrices g J instead of g: 

U 9 [f}(z)=U 9J [f}(z) = f((g T ) J z). (11) 

We claim that this new representation is in fact the (—p,j) representation under a change 
of variables. To see this, define a new set of homogeneous functions f(z) = f(Jz). The 
expression OU) for the Hermitian product is invariant under z — >■ Jz (the 2-form (TTUj) picks 
up a minus sign, but this is canceled as usual by a reversal of the integration interval). In 
terms of Jz, the transformation law (JTTT) takes the form: 

U g [f}(Jz) = f((g T yjz) = f(J(g T z)). (12) 

Therefore, the functions f(z) transform under the canonical rule (jSJ). Finally, the homogene- 
ity of f(z) is obtained from (CO) by interchanging A and A. The result is the homogeneity rule 
of the (p, —j) representation. We conclude that the action of g J on the (p,j) representation 
elements / is isomorphic to the action of g on the (p, —j) representation elements /. As is 
well known, the (p, —j) representation is isomorphic to the (—p,j) representation, so eq. (jUJ) 
follows. 

Coming back to the vertex amplitude ([3]), let us perform a change of integration variables 
at the nodes: 

9n^g J n - (13) 
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This does not affect the Haar measure, or the fixing of one of the g n to the identity. Therefore, 
dg n y dg n . Since ki = kf, the argument of the SX(2, C) character in ([3]) transforms as: 

%-(0%) ~* k i9s(i)(9t(i)) J '= (%«(»)%)) • (14) 
Using 06]), this implies: 

x^ihgs^g^)) ->■ x" 7JJ (%s(oV)) ■ ( 15 ) 

Since the entire operation was just a change of integration variables, we obtain the result: 

A v>J (hi) = A v ^(h t ) . (16) 

1. Alternative derivation without infinite- dimensional characters 

The above derivation of the symmetry (TIB"]) made free use of the infinite-dimensional 
SX(2,C) characters x 1 ^'^ ■ These are defined only in a distributional sense. Therefore, it's 
worth presenting an equivalent derivation which doesn't invoke these quantities. Let us 



consider the alternative definition of ([3]) in terms of representation matrix elements 131 ] : 



A*<i(h)= I n^IIE^' + 1 ) E Dj (h) m m 'D^(g s{l) g^y m ' jm , (17) 

where D^(h) is an SU(2) representation matrix in the magnetic-number basis, and D 1: >j(g) 
is an SL(2, C) representation matrix in the (spin, magnetic number) basis. Once again, we 
can perform the substitution g n — > g 3 n on the integration variables without changing the 
result. This substitution sends g s (i)g^ {.9s{i)9iq^) j ■ Then to demonstrate the symmetry 
( fl6l) . it suffices to establish the relation: 

D ^-(gV'-'j,,, = D^(g J y m ' jm , (18) 

for an arbitrary SL(2, C) element 9a(t)9fm = 9- To prove this relation, we will use the 
decomposition: 

g = Ul b(e)u 2 , (19) 
where Uj, u% G SU(2), and b(e) is a pure boost in the tz plane with boost parameter e: 

''(.)= r v i • <*» 



£ 



-1 
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The SU(2) elements U\,U2 are invariant under the J-conjugation. Also, the corresponding 
representation matrices are independent on 7, and are simply given by the spin-j SU(2) 
matrices: 

iy- j {«,)"'■[„„ = D-^( Ui y m ' jm = />'(»,)'"',„ . (21) 

It therefore suffices to demonstrate the relation ( ITS"]) for the pure boost b(e). The effect of 
J-conjugation on b(e) is to send the boost parameter to its inverse: e — > e . We must 
therefore show that: 

(b(e)y m ' jm = D™ {bie- 1 ))^' jm , (22) 

The matrix elements of a pure boost b(e) in the (p, j) representation are explicitly known 
jl^ . The elements of interest to us are obtained by setting j = j' = u in eq. (4.11) of [l^ . 
We get: 

D ii,3 (b(e)) jm ' jm = ^' £ 2(i+ m +i(i+n/2)) 

■f(i+j fl + ^V j + m + l; 2(j + l); 1 - ^ 



2 

where F(a, (3; x; y) is the hypergeometric function. The desired relation fl22|) can be derived 
directly from (123]) . using the following property of F(a, (3; x; y): 

F(a, (3; x; y) = (1 - y)^F (x - «, (3; x; ^) . (24) 

Reeling back the string of logic, we have thus established the symmetry ( fl6l) of the vertex 
amplitude. 



B. Derivation from semicoherent states 



In 12J], the vertex amplitude for a semicoherent state with spins j a b and spinor parameters 



£ a b is written as: 

Av(jab,U,Zba) = (-l) x / X\dg n FT Pab , (25) 

b=t(l) 

where n labels the graph's nodes, / labels the links, and a, b label the source and target nodes 
of each link. (— l) x is a global sign factor from combinatorics, and the "link propagator" 
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P ab is given by: 

Pab = P {g a 14>ab, 9b^4>ba) 

= c jt(2j b + 1) / - 2iPab t^iM ) 23ab (26) 

71 JCP 1 \\Z a b\\ 2 \\Zba\\ 2 \\\Z a b\\J \ \\Z a b\\\\Zba\\ J 



We refer the reader to [12| for a full explanation of this expression. The integral is performed 
over 2-spinors z a j, with the measure Q from ffTUj) . The Z spinors are a shorthand notation 
for Z a b = glz a b and Z ba = g\z a b- The angle brackets stand for the Hermitian inner product 
in C 2 , and the norms \\Z\\ are defined with respect to that inner product. The phase factor 
c ab is given by: 

jab + iPab , 

Cab = 1-2 T 2 • ( 2? ) 
VJab + Pab 

We set p a b = •yjab, which makes all the c a bS equal: 

Cab = c 7 = 7== • (28) 
V 1 + T 



Note the asymmetric treatment in ( 1261) of the link's source and target nodes a and b. We 

□ 

may interchange them by using the symmetry property of the (3 form [12]: 

/3(0 1; 2 ) = (-1)^/3(02,0!) . (29) 

We then have: 



= (-i) x / n (-i) 2 '^ (30) 

JSL(2,Q*-i n a=s(l) 
b=t{l) 

ha = P {9bZ4>ba, g a Z4>ab) 

c ab /f1 . , i\ f ^ f \\Zab\\\ 2lPab f ( — JZ a b,£,ab) (Zba,£,ba) Sx ~'' 



(2j a fe + 1) i n 7 i|2||y |2 \ II 7 
" JCP 1 H^aMI ll^all \|| z '6a||/ V ll^a&ll \\^ba\ 

Now consider the change of integration variables: 

g n g J n \ Zab Jzab ■ (32) 

The dg n and Q integrals are invariant under these changes; so is the fixing to the identity of 
one of the g n . The induced change in Z a b and Z ba is: 



J «6 



J Z a b', Zba — > J Zba ■ (33) 
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This turns (13TT) into: 



p v C ab /r) ■ , r\ f ^ /^ll^ll \ 2tPat f ( Z ab,£ab) {J Zba^ba)^^ 

^ ~{ Z]ab + ^ / 117 II2II7 112 My II 117 II II 7 II J ^> 

" -'CP 1 ll^a&ll II^Ml VII^Ml/ V H^afell ||^fea 



This is the same as P a b, except for the sign on JZ ba and an inversion of the (||Z ba ||/||Z afe ||) 2jp ' 
factor. This last difference corresponds to flipping the sign of p a b, or equivalently of 7. 
Keeping in mind also the dependence of c a b on p ab , we get: 

P fca ( 7 ) = (-l) 2 ^c 2 a6 P ab (- 7 ) • (35) 
The sign factor cancels nicely with the one in ( 13U1) . giving: 

where L is the number of links in the graph. There is a discrepancy between the symmetries 
( 1361 and ( TIB]) . It suggests that the vertex amplitudes fl3J and ( 1251) are not quite the same. 
Since ( f!6|) is a cleaner symmetry, it appears sensible to adopt a version of the semicoherent 
amplitude ( )25|) that respects it. This amounts to omitting the c a b factors in f[25]) - (|26|) and 
defining: 



A v {j ab ,U,£ba) = (-l) x / T\dg n TT P t 

JSL(2,C) N - 1 „ 7a 



SL(2,C)"-i V a"0 

&=t(0 (37) 

p l/ - , i\ f ^ /ll^&all f(Zab,£ab) (~ J Z ba , £,ba) \ 2j ^ 

r a b - -{^3ab + 1J / ||7 ||2|| „ ||2 I, y .I 1| y |||| 7 || 

JCP 1 Uracil ll Z M| \ll Z af>ll/ V II L ab\\ \\ ^ba\\ J 

This is the definition that we'll adopt in the numerical analysis below. As we'll see in 
section IIII D\ the omission of the c a b factors also leads to a more symmetric relation between 
parity-related amplitudes at fixed 7. 



C. Flipping 7 at each link separately 

The work reported in this section began with the hope that the vertex is not invariant 
under 7 — > —7. Then a modified vertex of the form A v>1 + A v _^ could have resolved the 
chirality problems. As we've seen, this is not the case, and the proposed modification is 
trivial. 

One may consider an alternative modification, which we now briefly discuss. 7 is tradi- 
tionally taken to be a global parameter. However, its actual usage in the vertex amplitude 
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(J3]) or (157)1 is to raise the representations on each link from SU(2) to SL(2, C) - an opera- 
tion which can be considered for each link separately. In particular, we can imagine using 
different signs of 7 on different links. To restore the symmetry between the links, we must 
then sum over all 2 L such possibilities. 

The arguments of subsections III Atill Bl which demonstrated invariance under a global 
flip of 7, no longer apply. Therefore, this "per-link" flipping of 7 may be a genuine and 
interesting modification of the model. However, we must note that in the large-spin limit 
studied in [12|, it degenerates back to the "global" flip A V}J + Ay_ 7 . Indeed, it was shown 
in 12J that to avoid exponential suppression, p a b must be proportional to jab with the same 
coefficient on all the links. Therefore, mixed terms with 7 on some links and —7 on the 
others won't contribute to the asymptotic amplitude. 

III. NUMERICAL RESULTS IN THE LARGE-SPIN 4-SIMPLEX LIMIT 
A. Introduction to the 4-simplex geometry 

We now turn to the parity structure of the large-spin 4-simplex amplitudes. We wrote 
a numerical script for the task. The script works with the semicoherent states described in 



121 ] , but without the c a b phase factors as we discussed in section III Bt The script is used to 
generate appropriate boundary-state parameters (j a b,^ab) from simpler geometric data, to 
find the two critical points (g a , z ao )± for the transition amplitude, to calculate the Hessian of 
the "action" at these points, and finally to calculate the "semicoherent" piece of the graviton 
propagator. 

In the Euclidean calculation of the graviton propagator the authors used a regular 4- 
simplex for maximal symmetry. In a Lorentzian signature, a regular 4-simplex doesn't exist. 
We therefore work with the next most symmetric possibility, i.e. an isosceles 4-simplex. It's 
composed of a regular "base" tetrahedron in the t = hyperplane centered at the origin, 
and 4 isosceles "side" tetrahedra whose apexes meet at a point (h, 0,0,0). In turn, each 
of the isosceles tetrahedra is composed of a regular base triangle (where they touch the 
base tetrahedron) and 3 isosceles side triangles (where they touch each other). The base 
tetrahedron is labeled as node no. 0, and the side tetrahedra are labeled as nodes 1 ... 4. 

Figure [U depicts (in one fewer dimension) the location of the tetrahedra in the "source" 3d 
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FIG. 1: A dimensionally-reduced depiction of the tetrahedra (here, triangles) that define the 4- 
simplex boundary state. The four (here, three) isosceles side tetrahedra are attached to the regular 
base tetrahedron from the inside. When the side tetrahedra are boosted by appropriate angles, 
the circled vertices join together. These boosts correspond to the SL(2, C) matrices g a at the two 
critical points. 

space, where the boundary state is defined. Here again, the base tetrahedron is centered at 
the origin. The side tetrahedra have a smaller height, and are glued to the base tetrahedron's 
faces from the inside (see section [III CI on this point). The SL(2,C) rotation matrices g a 
at the two critical points leave the base tetrahedron intact, while the side tetrahedra are 
boosted around their base triangles until their side triangles and apexes meet. The difference 
between the two critical points is that one boosts the side tetrahedra into positive t, and 
the other into negative t. 

We note that exchanging two tetrahedra among the equivalent set 1 ... 4 amounts to a 
parity-odd isometry on the 4-simplex. For instance, we can (and do) choose 4d axes so that 
tetrahedra 0,3,4 are centered on the z = hyperplane, while tetrahedra 1,2 are centered 
at opposite values of z. Then the exchange 1 -H- 2 corresponds to the parity-odd isometry 
z — > —z. This correspondence between parity and the exchange of two tetrahedra will play 
a crucial role in the analysis of the graviton propagator in section IIIIEI 

B. Description of the numerical calculation 

In this subsection, we describe the logical flow of the numerical calculation. Several 
details more related to programming than to geometry are deferred to the Appendix. 
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The recipe for semiclassical amplitudes in [12| begins with the parameters {j a b, Cab) of the 
boundary state. These are then translated into a discrete 4-simplex geometry by the critical 
points of the integrals (E7J), if the boundary parameters admit such critical points. The 
geometry at the critical points may be Euclidean, Lorentzian or degenerate (i.e. effectively 
3-dimensional). In the present context, our goal is slightly different. We wish to start 
with a boundary state which corresponds to a Lorentzian 4-simplex, and use this as a 
background for the graviton propagator. Therefore, the calculation works "backwards" : we 
start by specifying the desired geometry of the Lorentzian 4-simplex at the critical points, 
and calculate from that the required parameters {jabiCab) of the boundary state. 

The vertex positions of the base tetrahedron are hard-coded into the script. Its edge 
length is normalized to 1. This is the only hard-coded geometric data. All the rest is 
calculated dynamically in order to minimize human error, in particular with regard to signs 
and orientations. The script accepts as an input the rapidity with which the side tetrahedra 
are to be boosted at the critical points. From this we calculate their height and locations in 
the 3d source space, i.e. their vertex coordinates prior to boosting. 

Once the coordinates of the tetrahedra are known, the parameters {jab, Cab) of the bound- 



ary state are calculated, following the prescription of [12| in reverse (up to the caveat in sec- 
tion [TlLC])- We take the spins j a b to simply equal the corresponding triangle areas. The nor- 
malization is irrelevant for our purposes, and no harm is done by the non-half-integer values - 
in the stationary-phase limit, only the relative sizes of the spins matter. For a < b, we choose 
arbitrarily the phase of Cab according to the convention (e - ^/ 2 cos(0/2), e tLp ^ 2 sin(0/2)). The 
phase of Cba is then determined geometrically by the prescription in 12]. 

We calculate the rotation matrices at the critical points from the geometry of the 
4-simplex. The critical values z ab of the z a b spinors are then determined from the first 
critical-point equation in 121 ]. with the second component of each spinor normalized to 1. 

The major task, and the one which drove us to a numerical treatment in the first place, 
is to calculate the Hessian H of the "action" S [12] at the two critical points. For that 
purpose, we must first encode the integration variables {g a , z a b) into a set of 24 + 20 = 44 
non-redundant real quantities. Also, it's useful to choose these quantities so that the Haar 
and Q measures become trivial around the critical points. We do this as follows. For 
the rotation matrices, we first note that only 4 of them should be integrated over. We 
therefore keep g Q = 12x2, an d encode only the four others. For each of those, we encode 
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not g a itself, but the matrix g a = g a {gt) 1 '> which represents the deviation from the critical 
point. As the 6 independent real components of g a , we use the real and imaginary parts of 
(G?a)oi + (<7a)io)/2, i{{g a )oi ~ (<7 a )io)/2 and (<? a )oo- Near the identity, these correspond to 
the Cartesian components of the rotation and boost generators, which ensures even spacing 
under the Haar measure. To encode the z a b spinors, we simply use the real and imaginary 
parts of (z a b)o, after normalizing (z a b)i to 1. 

With these ingredients in place, we calculate the Hessian by straightforward numerical 
differentiation. We will use the determinant of H for the relative weights of the critical 
points in the amplitude, and its inverse for the graviton propagator. 

We next calculate the gradients of the metric insertion functions q^ b for a, b ^ n at the 
two critical points. The prescription for the metric insertions is taken from jg], with two 
trivial modifications: to accommodate for all the cases a < b, a = b and a > b, and to 
correct for the ab <-> ba discrepancy in conventions between jg] and 12]. This results in the 



following formula for q% (g a , z, 



<ii, 



^ab 



A, 



Una' < 



{ijnaf a = b 

A na • A n b a 7^ b 

((jZ nay Cna) 
(^noi Cna) 

(<j(JZ na ),£ na ) 



n < a 



(38) 



(39) 



n > a 



where a are the Pauli matrices. 

At each of the two critical points, we contract the gradients of q^ b with the inverse Hessian 
to obtain the semicoherent graviton propagator: 



Mab)(cd) _ / TT-\\ij f ab\l r cd\l 



(40) 



The propagator component Gmn represents the correlator between the metric elements 
ma -mb and nc-nd. The indices i,j in (j40p run over the independent components of (g a , z a b). 

Two caveats are in order. First, the script only handles the case m ^ n and a, b,c,d ^ 
m, n, when the double metric insertion is simply given by the product of two single insertions. 
Second, as mentioned in the Introduction, the quantity that we're computing is not the full 
graviton propagator, but only its "semicoherent" piece. The full propagator would have 
been obtained from a fully coherent boundary state, with a weighted sum over spins. This 
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sum would pick just one of the two critical points, and add a contribution to (|4"0|) from 
the derivatives with respect to j a b. In fact, this contribution is the one containing the 
standard Regge propagator. We ignore it in the present context, because we are interested 
in parity violation, and that comes from the derivatives with respect to (g a ,Zab)- Thus, we 
are calculating only the potentially troublesome part of the propagator. It can be cleanly 
separated from the j a b~ derivatives contribution, due to the vanishing of the Hessian elements 
H jg ,H jzj H jg [8|. 

We stress that all of the above can be done analytically. In particular, numerical dif- 
ferentiation is not necessary: not only the relevant functions, but also their derivatives are 
worked out analytically in [8]. The virtue of the numerics in this case is purely tactical - 
it provides a shortcut to some basic answers, allowing us to avoid tedious and error-prone 
derivations. Also, it will be able to provide a useful check for future analytical results. 
In fact, the hands-on numerics has already helped us uncover a minor problem in earlier 
analytical work, as we now describe. 



C. A note on the orientation of the side tetrahedra 



The critical-point equations on (g a , z a ^) do not have solutions for generic boundary states 
(jab,£,ab)- In order to correspond to a semiclassical Regge geometry, the boundary-state pa- 
rameters must satisfy certain constraints. These constraints have a geometric interpretation, 
as described respectively in ljj and 12[ for the Euclidean and Lorentzian models. In the 
Euclidean, the £ {, should be understood as the spinor "square roots" of the outgoing face 
normals of the tetrahedra, embedded in 3d space. Then the critical-point solutions corre- 
spond to folding these tetrahedra into a 4-simplex jusing 4d rotations. The same procedure 
was then carried over to the Lorentzian model in jl2[. As we noticed while looking for crit- 
ical points in the numerics, there is in fact a subtle difference between the two cases which 
must be taken into account. 

For concreteness, consider the isosceles 4-simplex described in sections IHI AtlllTBl In the 
Euclidean model, its tetrahedra can be embedded in 3d space as a four-pointed "star" : the 
base tetrahedron is situated at the center, with the side tetrahedra glued to its faces from 
the outside. To fold this into a 4-simplex, the side tetrahedra must be rotated along their 
heights into the 4th dimension by some obtuse angle. Now, in the Lorentzian case such 



15 



a rotation is not continuous with the identity: it would require passing through null and 
timelike configurations. In other words, it isn't part of SL(2, C). To fix this, we must glue 
the side tetrahedra to the inside of the base tetrahedron, as in figure [TJ That way, they can 
be folded together by a finite boost. 



Now, the critical-point equations in [12| require the spinors and J^ a to point in 
the same direction after the 4d rotations. In other words, £ a & and £& a should point in 
opposite directions. This is the case in the Euclidean setup, with all the £'s pointing along 
the outgoing normals of their respective tetrahedra. However, now that we took the side 
tetrahedra to lie inside the base tetrahedron, the base and side outgoing normals are parallel 
rather than antiparallel. This can be cured by taking the £'s for the side tetrahedra (or the 
base tetrahedron) to point along the ingoing normals to their respective faces. 

This conclusion can be generalized to arbitrary Lorentzian 4-polytopes with spacelike 
boundaries. The boundary polyhedra can be classified according to the time-orientation 
of their outgoing 4d normal (not to be confused with the face normals discussed above). 
In constructing the semiclassical boundary state, we must take the £ spinors for the past- 
pointing polyhedra (say) to point along the outgoing face normals, and for the future- 
pointing polyhedra - along the ingoing face normals. This was the prescription we used in 
the numerical script. 

D. The coefficients of the two critical points 



In 12j . the asymptotic analysis of the Lorentzian 4-simplex amplitude (with Lorentzian 



critical points) ends with the following expression: 

A v = (-1)*+ M (N+e iS * e °° a + N-e- 43 **"') , (41) 
where the sign factor comes from combinatorics, and SR egge is the Regge action: 

SRegge = 7 3ab®ab , (42) 
a<b 

with Q ab the dihedral angles (actually, boost parameters). The coefficients N± of the two 
parity-related critical points are given by: 
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where the piece in the second parentheses is the one which depends on the critical point. 
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We again omitted the c a b phase factors, as discussed in section III Bl The authors of 
left the ratio N + /N_ as an open question. We studied this question using the numerical 
script. Specifically, the script calculates the quantity (^JdetH Y[ \\Z a b\\ 2 \\Z ba \\ 2 ^j at the 
two critical points, keeping in mind that the Q measure is trivial in our conventions. The 
numbers point clearly towards the simple result: 

N + = NZ . (44) 

In other words, the amplitude (T4T|) is real. In particular, the magnitudes \N+\ = |AL| are 
equal. Furthermore, at small 7 we have: 

7 ^0 N + = N^=N^=JT. (45) 

Here and below, we present extrapolations from numerical results as exact analytical re- 
lations. In all such cases, what's actually implied is that the LHS and RHS agree to an 
accuracy of several digits, for a representative sample of points in the parameter space. 

The results ()44l) -( l45l) may be compared with the Euclidean situation [ll], where the 
coefficients of exp(iSR egg e) and exp(—iSR egge ) are equal for arbitrary 7 (and two additional 
terms appear). Note that the simple relation ( jUJ) is sensitive to the overall phase of the 
boundary state. It holds when the overall phase of and £& Q is fixed by the 3d rotation 
which makes the corresponding triangles congruent, as prescribed in [l2j. The result is also 
dependent on our omission of the c a b phase factors in (1371) and (1431) . 

As a cross-check, the numerical results for N± confirm the 7 — > —7 symmetry claimed in 
section [TTl Concretely, the identity N±(j) = N T (— 7) is seen to hold. 



E. The semicoherent graviton propagator 

Our main object of interest is the graviton propagator, in particular its parity-related 
components. More specifically, as discussed in section IIIIBt we are calculating the "semi- 
coherent" propagator fT40l) . where the summed indices i,j run over the independent com- 
ponents of (g a , z a b) and not over the spins j a b- We calculate the quantity ( 140]) at the two 
critical points separately, noting that only one will remain after the sum over spins. Also, 
we reiterate that ( 1401) is valid only for m 7^ n and a, b, c, d 7^ m, n. 
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As a first result of this calculation, we again see a confirmation of the 7 — > —7 symmetry, 
i.e. G^ )M (7) = G^ )M (-7)- To go further, we restrict our attention to the components 
G^ ){cd \ with 1 < n, a,b } c } d < 4 and a, b,c } d 7^ n. As expected from rotational invariance, 
all even permutations of the equivalent nodes (1, 2, 3, 4) result in the same value for G^® . 
However, odd permutations, which are related to the even ones by spatial parity, result in 
a different value. As representatives of these two equivalence classes, we pick G^f 1 ^ and 
Gg 4)(23) . We now evaluate these components at the same critical point (the "positive" one), 
since that's how they will enter the full coherent propagator. In figure |2j we plot the real 
and imaginary parts of these quantities (divided by 7 4 ), as well as their complex phase, as 
functions of 7. We use fixed values for the spins and a fixed boost parameter determining 
the shape of the 4-simplex. From the numbers and from the graphs, we make the following 
observations: 

1. The ratio of the two parity- related propagator components is: 

^i 3)(24) /4i 4)(23) =e 2m/3 , (46) 

exactly as in the Euclidean case Q. In particular, we see that this result holds also 
for 4-simplices which are not regular. 

2. At small 7, the semicoherent propagator scales as 7 4 . This can be seen from the fact 
that G/7 4 remains regular as 7 — > 0. The same scaling law has been found to hold in 
the Euclidean model [7j. Since the P-invariant Regge piece of the propagator scales 

it will dominate at small 7. 



as 7 3 
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3. At small 7, the two parity-related components become each other's conjugates: 



7^0 =► Gg 3){24) = Gg 4)(23) . (47) 

The same was found in [3] for the Euclidean model. 

PI 

4. In the Euclidean regular 4-simplex [71], the semicoherent piece of the propagator is 



proportional to 7 4 (7 — iy/TEj). With such a dependence, the graphs in figure 2(a) 
would have been straight lines. We conclude that the Lorentzian propagator does not 
have the same polynomial dependence on 7. 
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0.003 




7 



(a)The real and imaginary parts of two parity-related propagator components. We divide by 7 4 , 
because the corresponding graphs for the Euclidean regular 4-simplex are straight lines. 




0.0 0.5 1.0 1.5 2.0 2.5 3.0 

7 

(b)The complex phase of the same two propagator components. 

FIG. 2: The "semicoherent" piece of the 4-simplex graviton propagator at one of the critical points 
as a function of the Immirzi parameter. The boost parameter of the side tetrahedra is fixed to 1. 
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Finally, another regularity emerges from comparing the propagator components at the two 
critical points. Specifically, a propagator component at one point equals the complex con- 
jugate of the parity-related component at the other point. For example: 



G (23)(24)\ = ( G m(23)\ (4g) 



; + - r oi ; 

We stress that the above results for the graviton propagator do not depend on the choice 
of phases for the spinors defining the boundary state. This is in contrast to the results 
for the transition amplitude in section IIII Dl The reason for this robustness is that the 
propagator is a ratio [7|, |8j between two amplitudes - one with the metric insertions and one 
without. Any phase factors in the boundary state cancel when taking this ratio. Thus, it 
appears unlikely that the parity-violating result fj4*6]) is an artifact of any phase or orientation 
choices in the definition of the boundary state (other than the global orientation, which is 
the subject of interest). 



IV. DISCUSSION 



We have seen that the Lorentzian 4-simplex graviton propagator is not parity- invariant, 
exactly as in the Euclidean case. We've also shown that flipping the sign of 7 has no effect 
on the spinfoam vertex formula, so it cannot help to resolve the problem. If LQG is to be 
taken seriously as reproducing GR at large distances, this issue must be addressed. 

Perhaps the problem will go away once larger graphs and spinfoams are taken into ac- 
count. Such calculations should be performed in order to find out. Statistically, one may 
expect large graphs with small spins on the links to dominate over small graphs with large 
spins. Therefore, a good place to start may be the single-vertex approximation for a large 
graph with spin- 1/2 links. Perhaps the spin-1/2 4-simplex will be a good warmup exercise, 
though it isn't directly relevant to the large-distance limit. 

Another possibility is that the correlator Gmn between two single nodes is not the 
correct quantity to describe the large- distance graviton propagator. In the crude 4-simplex 
approximation, there's really no other choice. However, given a large boundary graph with 
many Planck-scale nodes, we may consider a correlator that's smeared over many nearby 
nodes around each of the two endpoints. In fact, one might say that calculating a large- 
distance propagator between points defined at a Planckian resolution is outright suspect. 
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We should note here that the parity problem, in both the Euclidean and Lorentzian 
models, goes away in the limit 7—7-0. In fact, the evidence so far seems to indicate 



that the large-spin, small-7 limit precisely reproduces classical Regge gravity [8|, ll5|, [16|. 
We are not sure what to conclude from this. One difficulty should certainly be noted: 
in the aforementioned limit, Regge gravity is indeed reproduced precisely, e.g. with no 
renormalization of Newton's constant. This stands in tension with the LQG black-hole 



entropy calculations 17H19|. where agreement with the Bekenstein formula demands either 
an order-1 value for 7 or a 7-proportional renormalization ratio Gm/Guv- 

It is possible that none of the above solves the issue, and that the theory must be modified. 
As we've seen, the vertex formula already appears P-invariant. Therefore, the necessary 
modification may lie in the geometric interpretation of the boundary states. Rescaling the 
graviton field components to obtain a P-symmetric propagator falls within this category. It 
is an extreme possibility, as it amounts to changing the time-honored kinematics of LQG. 
Finally, it may be that the problem will be solved by large graphs with small spins together 
with the per-link modification proposed in section III CI 

Full analytical understanding of the 4-simplex system is an important task in itself. 
Our numerical observations P4|) - fj45|) and fj46|) - fj48|) can and should be derived analytically. 
At present this seems like tedious work. Given the simplicity of the final results, we are 
probably missing some crucial ideas. As an intermediate step, symbolic math software may 
be considered. 

The numerical 4-simplex script presented here can be used for questions other than parity 
invariance. More generally, LQG is bursting with fundamental questions for which numer- 
ical answers will be vastly better than nothing. We should look forward towards more 
comprehensive numerical tools, which are still largely missing. 
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Appendix A: Technical comments regarding the numerical script 

1. Methods of invoking the script 

In writing the numerical script, efficiency was consistently sacrificed for simplic- 
ity. The 4-simplex calculations, with textual output, are carried out by the script file 
spinf oam_4s implex . py. The plots in figure |2]were produced using the script file G_plot . py. 
The main script spinf oam_4s implex . py can be run as a standalone (see the usage_string() 
function), or from another program using the calc_4simplex() function. A 'verbose' bi- 
nary flag causes the script to print out in detail the parameters of the boundary state and 
the critical points. A 'debug' flag prints out some consistency checks: that the critical- 
point equations are satisfied, that the gradient of S at the critical points is small, etc. A 
'crit-points' option can instruct the script to process just one critical point, as is relevant 
for the fully coherent boundary state, or none, if we are only interested in the boundary 
state parameters. 

2. Vertex placement and symmetry planes of the base tetrahedron 

In section IIII B\ we mentioned that the base tetrahedron's vertices are hard-coded so as 
to obtain a regular tetrahedron with edge length 1 centered at the origin. This leaves us 
with the freedom to choose the orientation of the tetrahedron. Our choice was to place one 
of the face normals along the x axis, and to align one of the symmetry planes with the xy 
plane. In making this choice, it was important that none of the face normals are parallel to 
the z axis - this is the only failsafe against singularities due to a vanishing component of a 
2-spinor. 

3. Choice of differentiation step 

In our calculation of the gradients and the Hessian at the critical points, the differentiation 
step is hard-coded. The optimal order of magnitude for the step was found by a trial-and- 
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error process. There are two opposing constraints at work here. On one hand, the step 
should be large enough to be well above the floating-point accuracy. For this, we ensure 
that the results of differentiation are insensitive to order-of-magnitude changes in the step. 
On the other hand, the step should be small enough to avoid smearing the results. For this, 
we ensure the smallness of the "action" gradients at the critical points. We found a range 
of several orders of magnitude which satisfies both requirements, and chose a step near the 
middle of that range. 

This manual choice of the step should be kept in mind (and hopefully improved upon) if 
one wishes to adapt the code for a different calculation. 
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